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Wo present the first measurement of the integrated forward-backward cliargc asymmetry in top- 
antitop quark pair (tt) production in proton-antiproton (pp) collisions in the lepton+jets final state. 
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Using a 6-jet tagging algorithm and kinematic reconstruction assuming tt + X production and decay, 
a sample of 0.9 fb~^ of data, collected by the DO experiment at the Fermilab Tevatron Collider, is 
used to measure the asymmetry for different jet multiplicities. The result is also used to set upper 
limits on tt + X production via a Z' resonance. 

PACS numbers: 12.38.Qk, 12.60.-i, 13.85.-t, 13.87.Ce 



At lowest order in quantum chromodynamics (QCD), 
the standard model (SM) predicts that the kinematic 
distributions in pp — > ti + X production are charge 
symmetric. But this symmetry is accidental, as the 
initial pp state is not an eigenstate of charge conjuga- 
tion. Next-to- leading order (NLO) calculations predict 
forward-backward asymmetries of (5-10)% [11, 0|, but re- 
cent next-to-next-to-leading order (NNLO) calculations 
predict significant corrections for tt production in associ- 
ation with a jet The asymmetry arises mainly from 
interference between contributions symmetric and anti- 
symmetric under the exchange t t [I], and depends on 
the region of phase space being probed and, in particu- 
lar, on the production of an additional jet 2]. The small 
asymmetries expected in the SM make this a sensitive 
probe for new physics 4]. 

A charge asymmetry in pp — > tt + X can be ob- 
served as a forward-backward production asymmetry. 
The signed difference between the rapidities [51 of the 
t and i, Ay = yt ~ yt, reflects the asymmetry in tt pro- 
duction. We define the integrated charge asymmetry as 
Afb = {N{- Nh) / {Nf + Nh), where A^'f (A"b) is the num- 
ber of events with a positive (negative) Ay. 

This Letter describes the first measurement of A[h in 
pp tt + X production. The 0.9 fb~^ data sample used 
was collected at ^/s ~ 1.96 TeV with the DO detector 0, 
using triggers that required a jet and an electron or muon. 
In the lepton-|-jets final state of the ti system, one of the 
two W bosons from the tt pair decays into hadronic jets 
and the other into leptons, yielding a signature of two b- 
jets, two light-flavor jets, an isolated lepton, and missing 
transverse energy {fr)- This decay mode is well suited 
for this measurement, as it combines a large branching 
fraction (~ 34%) with high signal purity, the latter a con- 
sequence of requiring an isolated electron or muon with 
large transverse momentum (px)- The main background 
is from VF-l-jets and multijet production. This channel 
allows accurate reconstruction of the t and t directions 
in the collision rest frame, and the charge of the electron 
or muon distinguishes between the t and t quarks. 

The dependence of Afb on the region of phase space, as 
calculated by the mc@NLO event generator Q, is demon- 
strated in Fig. [T] The large dependence on the fourth- 
highest jet Pt is not available in the calculations of Refs. 
PiSSIj as these do not consider decays of the top quarks, 
and include only acceptance for jets from additional ra- 
diation. 

We conclude that acceptance can strongly affect the 
asymmetry. To facilitate comparison with theory, the 
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analysis is therefore designed to have an acceptance 
which can be described simply. Event selection is lim- 
ited to either: (i) selections on directions and momenta 
that can be described at the particle level (which refers 
to produced particles before they start interacting with 
material in the detector) or (ii) criteria with high signal 
efficiency, so that their impact on the region of accep- 
tance is negligible. In addition, the observable quantity 
and the fitting procedure are chosen to ensure that all 
events have the same weight in determining the asym- 
metry. 

The measurement is not corrected for acceptance and 
reconstruction effects, but a prescription provides the ac- 
ceptance at the particle level. Reconstruction effects are 
also accommodated at the particle level by defining the 
asymmetry as a function of the generated | Ay \ : 



^fb(|Ay|) 



g{\Ay\) ^ gHAy\) 
g{\Ay\)+g{-\Ay\)' 



(1) 



where g is the probability density for Ay within the ac- 
ceptance. This asymmetry can be folded with the "geo- 
metric dilution," T), which is described later: 



^fb 



red 



Atb (Ay) V {Ay) [g {Ay) + g {-Ay)] dAy. 

(2) 

This procedure yields the predictions in Table [J The 
values are smaller than those of Ref. P, 0], because of 
the inclusion of jet acceptance and dilution. 

We select events with at least four jets reconstructed 
using a cone algorithm [sl] with an angular radius TZ — Q.h 
(in rapidity and azimuthal angle). All jets must have 
Pt > 20 GeV and pseudorapidity (relative to the re- 
constructed primary vertex) \r]\ < 2.5. The leading jet 
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TABLE I: Predictions based on mc@nlo. 







AZ" (in %) 


> 4 
4 

^ 5 


0.8±0.2(stat. 
2.3±0.2(stat. 
-4.9±0.4(stat. 


) ±1.0(accept.) ±0.0(dilution) 
) ±1.0(accept.) ± O.l(dilution) 
) ±1.0(accept.) ± 0.2(dilution) 



TABLE IL Parameters of the dilution. The ±1 sd values 
include both statistical and systematic uncertainties. 



must have pt > 35 GeV. Events are required to have 
> 15 GeV and exactly one isolated electron with 
Pt > 15 GeV and \ri\ < 1.1 or one isolated muon with 
Pt > 18 GeV and |?7| < 2.0. More details on lepton iden- 
tification and trigger requirements are given in Ref. 
Events in which the lepton momentum is mismeasured 
are suppressed by requiring that the direction of the ^Jt 
not be along or opposite the azimuth of the lepton. To 
enhance the signal, at least one of the jets is required 
to be identified as originating from a long-lived b hadron 
by a neural network 6-jet tagging algorithm The 
variables used to identify such jets rely on the presence 
and characteristics of a secondary vertex and tracks with 
high impact parameter inside the jet. 

The top quark pair is reconstructed using a kinematic 
fitter [11], which varies the four- momenta of the de- 
tected objects within their resolutions and minimizes a 
statistic, constraining both W boson masses to ex- 
actly 80.4 GeV and top quark masses to exactly 170 GeV. 
The 6-tagged jet of highest pT and the three remaining 
jets with highest px are used in the fit. The 6-tagging 
information is used to reduce the number of jet-parton 
assignments considered in the fit. Only events in which 
the kinematic fit converges are used, and for each event 
only the reconstruction with the lowest is retained. 

The jet-pT selection criteria strongly affect the ob- 
served asymmetry (see Fig. [T]), and this must be con- 
sidered when comparing a model to data. Fortunately, 
these effects can be approximated by simple cuts on 
particle-level momenta without changing the asymmetry 
by more than 2% (absolute). This is verified using sev- 
eral simulated samples with generated asymmetries and 



particle jets clustered using the pxcone algorithm [12 1 
("E" scheme and TZ — 0.5). The particle jet cuts are 
Pt > 21 GeV and |r/| < 2.5, with the additional require- 
ment on the leading particle jet pT > 35 GeV and the 
lepton requirements detailed above. Systematic uncer- 
tainties on jet energy calibration introduce possible shifts 
of the particle jet thresholds. The shifts are GeV for 
the leading jet and ^\'^GeV for the other jets, for ±1 
standard deviation (sd) changes in the jet energy calibra- 
tion. The resulting changes in the asymmetry predicted 
using McOnlo are of the order of 0.5%. The eflFect of 
all other selections on the asymmetry is negligible. The 
predictions in Table [J use a more complete description of 
the acceptance based on efficiencies factorized in pT and 
r], accurate to < 1% (absolute). 

Misreconstructing the sign of Ay dilutes the asymme- 
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FIG. 2: The geometric dilution and its uncertainty band 
as a function of generated \Ay\ for standard model ti + X 
production and > 4 jets. 



try. Such dilution can arise from misidentifying lepton 
charge or from misreconstructing event geometry. The 
rate for misidentification of lepton charge is taken from 
the signal simulation and verified using data. False pro- 
duction asymmetries arising from asymmetries in the rate 
for misidentification of lepton charge are negligible ow- 
ing to the frequent reversal of the DO solenoid and toroid 
polarities. 

The dilution, V, depends mainly on \Ay\. It is defined 
as P = 2P — 1 , where P is the probability of reconstruct- 
ing the correct sign of Ay. It is obtained from tt + X 
events generated with pythia [13| and passed through a 
GEANT-based simulation [l3| of the DO detector, and is 
parametrized as: 



V i\Ay\) = Co In 1 + ci |Ay| + C2 \Ay[ 



(3) 



with the parameters given in Table HIl (see Fig. [5]). 

As this measurement is integrated in \Ay\, the depen- 
dence of the dilution on \Ay\ introduces a model de- 
pendence into any correction from observed asymmetry 
{A'^^) to a particle-level asymmetry. Such a correction 
factor would depend not only on the model's \Ay\ distri- 
bution, but also on its prediction of Af^ {\Ay\). Further- 
more, such a correction would be sensitive to small new 
physics components of the selected sample. We therefore 
present a measurement uncorrected for reconstruction ef- 
fects and provide the reader with a parametrization of V 
that describes these effects, to be applied to any model. 
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The dilution depends weakly on other variables corre- 
lated with Afb, such as the number of jets. This possible 
bias is included in the systematic uncertainties. Non- 
standard production mechanisms can affect reconstruc- 
tion quality, primarily due to changes in the momenta of 
the top quarks. By studying extreme cases, we find that 
when comparing non-standard tt + X production to data 
an additional 15% relative uncertainty on Af^ is needed. 

The main background is from VF-|-jets production. To 
estimate it, we define a likelihood discriminant C using 
variables that are well-described in our simulation, pro- 
vide separation between signal and W^-l-jets background, 
and do not bias \Ay\ for the selected signal. The follow- 
ing variables are used: the pt of the leading ^-tagged jet, 
the statistic from the kinematic fit, the invariant mass 
of the jets assigned to the hadronic W boson decay, and 
k^^^ = p^i'^i?™*", where it!'"'" is the smallest angular dis- 
tance between any two jets used in the kinematic fit, and 
p™'" is the smaller of the corresponding jets' transverse 
momenta. 

The next largest background after VF+jets is from 
multijet production, where a jet mimics an isolated elec- 
tron or muon. Following the procedure described in Ref. 
0, the distributions in likelihood discriminant and re- 
constructed asymmetry for this background are derived 
from samples of data that fail lepton identification. The 
normalization of this background is estimated from the 
size of those samples and the large difference in efficien- 
cies of lepton identification for true and false leptons. 
The effects of additional background sources not consid- 
ered explicitly in extracting A^^; namely Z-|-jets, single 
top quark, and diboson production; are evaluated using 
ensembles of simulated datasets and found negligible. 

The sample composition and Afb are extracted from a 
simultaneous maximum-likelihood fit to data of a sum 
of contributions to C and to the sign of the recon- 
structed Ay (Aj/roco) from forward signal, backward sig- 
nal, VF-l-jets, and multijet production. Both signal con- 
tributions are generated with PYTHIA, have the same dis- 
tribution in C, and differ only in their being reconstructed 
as either forward or backward. The T4^-|-jets contribution 
is generated with ALPGEN [l^ interfaced to pythia and 
has its own reconstructed asymmetry. Although W bo- 
son production is inherently asymmetric, the kinematic 
reconstruction to the tt + X hypothesis reduces its recon- 
structed asymmetry to [4.4 ± 1.6 (stat.)] %. The multijet 
contribution is derived from data, as described above. 
The fitted parameters are shown in Table IIIII Correla- 
tions between the asymmetry and the other parameters 
are < 10%. The fitted asymmetries in data are consistent 
with the SM predictions given in Tabic [D In Fig. [3] we 
compare the fitted distributions to data for events with 
> 4 jets. 

The dominant sources of systematic uncertainty for the 
measured asymmetry are the relative jet energy calibra- 
tion between data and simulation (±0.5%), the asymme- 



TABLE III: Number of selected events and fit results in data. 





^ 4 Jets 


4 Jets 


^ 5 Jets 


No. Events 


376 


308 


68 


tt + X 

W-l-jets 

Multijets 


2661^^ 
70±21 
40±4 

(12±8)% 


214±20 
51+19 

"^-18 

32.7133 

(19±9)% 


54li^ 
7+11 
' -5 

Y 1+1-6 

(-16+!?)% 



try reconstructed in W^-|-jets events (±0.4%), and the 
modeling of additional interactions during a single pp 
bunch crossing (±0.4%). The total systematic uncer- 
tainty for the asymmetry is ±1%, which is negligible 
compared to the statistical uncertainty. 

We check the simulation of the production asym- 
metry, and of the asymmetry reconstructed under the 
tt + X hypothesis in the M^+jets background, by re- 
peating the analysis in a sample enriched in VF+jets 
events. The selection criteria for this sample are iden- 
tical to the main analysis, except that we veto on any 
6-tags. Both the fully reconstructed asymmetry and the 
forward-backward lepton asymmetry are consistent with 
expectations. We also find that the fitted sample com- 
position (Table IIII|) is consistent with the cross section 
for tt + X production obtained in a dedicated analysis 
on this dataset. We check the validity of the fitting pro- 
cedure, its calibration, and its statistical uncertainties 
using ensembles of simulated datasets. 

To demonstrate the measurement's sensitivity to new 
physics, we examine tt production via neutral gauge 
bosons {Z') that are heavy enough to decay to on-shell 
top and antitop quarks. Direct searches have placed lim- 
its on tt production via a heavy narrow resonance [13], 
while the asymmetry in tt production may be sensitive 
to production via both narrow and wide resonances. The 
Z' tt channel is of interest in models with a "leptopho- 
bic" Z' that decays dominantly to quarks. We study the 
scenario where the coupling between the Z' boson and 
quarks is proportional to that between the Z boson and 
quarks, and interference effects with SM tt production 
are negligible. Using pythia we simulate tt production 
via Z' resonances with decay rates chosen to yield nar- 



row resonances as in Ref. jl7l |. and find large positive 



asymmetries [(13-35)%], which are a consequence of the 
predominantly left-handed decays. We predict the distri- 
bution of Afb as a function of the fraction (/) of tt events 
produced via a Z' resonance of a particular mass from 
ensembles of simulated datasets. We use the procedure 
of Ref. [18] to arrive at the limits shown in Fig. [H These 
limits can be applied to wide Z' resonances by averaging 
over the distribution of Z' mass. 

In summary, we present the first measurement of the 
integrated forward-backward charge asymmetry in ti+X 
production. We find that acceptance affects the asym- 
metry and must be specified as above, and that correc- 
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FIG. 3; Comparison of data for > 4 jets with the fitted model as a function of £. for events reconstructed (a) as forward 
(Aj/reco > 0) and (b) as backward (Ayreco < 0). The number of events from each source is listed with its statistical uncertainty. 
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FIG. 4: 95% C.L. limits on the fraction of tt produced via a 
Z' resonance as a function of the Z' mass, under assumptions 
detailed in the text. Limits expected in the absence of a Z' 
resonance are shown by the dashed curve, with the shaded 
bands showing limits one and two standard deviations away. 
The observed limits are shown by the solid curve, and the 
excluded region is hatched. 



tions for reconstruction effects are too model-dependent 
to be of use. We observe an uncorrected asymme- 
try of A'^" = [12 ± 8 (stat.) ± 1 (syst.)] % for tt + X 
events with ^ 4 jets that are within our acceptance, 
and we provide a dilution function (Eq. [J]) that can 
be applied to any model (through Eq. For events 
with only four jets and for those with ^ 5 jets, we 
find A°i^^ = [19 ± 9 (stat.) ± 2 (syst.)] % and A"^^' = 
[— 16l}y (stat.) ± 3 (syst.)] %, respectively, where most 
of the systematic uncertainty is from migrations of events 
between the two subsamples. The measured asymmetries 
are consistent with the mc@nlo predictions for standard 
model production. 
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